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Abstract. In computer algebra there are different ways of approaching the 
mathematical concept of functions, one of which is by defining them as so- 
lutions of differential equations. We compare different such approaches and 
discuss the occurring problems. The main focus is on the question of deter- 
mining possible branch cuts. We explore the extent to which the treatment of 
branch cuts can be rendered (more) algorithmic, by adapting Kalian's rules to 
the differential equation setting. 
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1. Introduction 

In Mathematics, the standard definition of a function is given, e.g., by Bour- 
baki [Bou68, §3.4]: it consists of the domain and codomain, and requires that a 
function be total and single-valued. Mathematical practice is usually looser, and 
tends to define functions locally, or "in a suitable open subset of C" . Usually, there- 
fore, users of computer algebra systems will use notations such as "ln(x)" without 
bothering too much about the actual domain and codomain of the function In they 
have in mind, hoping that the designers of the system have implemented what they 
need. An essential and important problem in any such definition is the determina- 
tion of possible branch cuts (also known as slits), see [DF94]. 

A simple way to define a function is by composing previously defined functions, 
provided their domains and codomains are compatible. This is perhaps best exem- 
plified by 

'1 



\fz := exp ^— In z 

which implies that the branch cuts of s/z are inherited from the definition of In z. 

The challenge is, of course, to decide which formula we should take. For instance, 
as reported in [Kah87, pp. 210-211], the classic handbook [AS64] changed its 
interpretation of the branch cuts for arccot from the original first printing (here 
denoted as arccoti) to that of the ninth (and subsequent) printings, and [fSTIO] 
(here denoted as arccotg). Both versions of arccot are related to arctan by the 
formulae: 

arccoti (x) = 7r/2 — arctan(aj), 

arccotg (x) = arctan(l/a;). 

The rationale for the "correct" choice is discussed further in Section 3.2. 

Difficulties occur also when dealing with functions without having precise def- 
initions for them. One of the early successes of computer algebra was symbolic 
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integration. This area is largely based on differential algebra [Bro05]. There, the 
expression ln(a;) is not even a function, but some element 9 in a differential field, 
such that 9' = 1/x. But both ln(x) and — ln(l/x) have the same derivative so 
that their difference is a constant in the sense of differential algebra: its derivative 
is 0. However, with the definition of In from [AS64], the function ln(x) + ln(l/x) 
is for all x in C except for real negative x, where it is 2iir. Thus early versions 
of integrators would wrongly compute J_ 2 2xdx/(x 2 — 1) as by first computing 
correctly the indefinite integral as ln(x 2 — 1) and then subtracting its values at both 
end points. 

The definition of functions is obviously relevant to identities between functions: 
the functions involved should have, as a minimum, the same domain and codomain. 
For instance [Kah87, pp. 187-188], the function 

(1) g(z) :— 2 arccosh ^1 + — arccosh 

is not the same as the ostensibly more efficient 



(2) q(z) := 2 arccosh ^2(z + 3) 

unless they are defined over an identical domain that avoids the negative real axis 
and the area 

[z = x + iy:\y\< ^ X + ^ A -9/2 < x < - 3 j 

(and an identical codomain). Clearly, this last statement itself depends on the 
definitions taken for the functions arccosh and ^J. Here, the function is defined 
over C \ R - , while arccosh is defined over C \ (1 + R~). In other words, these 
functions have branch cuts located on horizontal left half-lines. 

In general, the location of branch cuts for classical functions follows very much 
from the mathematical tradition, as recorded in tables like [AS64]. We discuss 
here the possibility of automating the choice of these locations for a large class 
of functions "defined" by linear differential equations. We show that this is im- 
possible in general. Therefore, we consider a heuristic approach that gives correct 
results for the special functions listed in [AS64] and is applicable to newly encoun- 
tered functions. As a guide, we use Kahan's discussion of this question for inverse 
trigonometric functions [Kah87]. 

We are concerned here with automating the choices made, historically by table- 
makers, and now by the compilers of resources like [fSTIO] or the authors of systems 
such as [DDMF] , for functions defined by linear differential equations and their inital 
conditions. A particular context for the use of the functions may impose special 
constraints, and may even require different choices at different points of the same 
application [DavlO, section 4(ii)], but that is a different issue. 

2. Analytic Continuation 

In complex analysis, analytic functions are often defined first on a small domain 
and then extended by analytic continuation. Two large families of examples are 
discussed below: inverse functions and solutions of linear differential equations. 
The basic property underlying this approach is that two analytic functions defined 
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over a connected domain and coinciding over an open subset of it coincide over 
the whole domain and thus are identical (here we assume the codomain to be C). 
Another way of discussing branch cuts is thus in terms of connected domains where 
the function of interest is to be defined. 

2.1. Riemann Surfaces. A radical approach is to use Riemann surfaces as do- 
mains. They are maximal connected surfaces where the function is analytic. While 
theoretically appealing, the use of Riemann surfaces is not trivial in a computer 
algebra context (see, e.g., [Hoe05]). Since the domain is not a subset of C, but 
paths in the complex place, an ad hoc language for specifying its elements has to 
be designed. One possibility is to restrict the domain to paths that are piecewise- 
straight lines starting from the origin. This is also the approach taken in the 
Dynamic Dictionary of Mathematical Functions [DDMF] . Such a path is specified 
as a list of its "vertices" , for example, (0, 1 + i, 2, 1 — i, 0) denotes a diamond-shaped, 
clockwise path around 1. Thus for instance, one can define so that it takes the 
value 1 at (0,1), while it is equal to —1 at (0, 1, i, — 1, — i, 1) and to 1 again at 
(0,l,i,-l,-i,l,i,-l,-i,l). 

2.2. Positioning. In many applications however, users are interested in restricting 
the domain of their functions to the complex plane or a subset of it. In that 
case, the role of branch cuts is to define a connected domain where the function is 
analytic. Where to put the branch cuts is the positioning question. Apart from the 
connectivity and analyticity constraints and as long as only one function is involved, 
the location of the branch cuts is quite arbitrary. The situation is completely 
different as soon as several functions are involved and identities arc considered: the 
domains have to coincide. Thus branch cuts have to be chosen in a consistent way 
inside a corpus of functions of interest. 

2.3. Adherence. It is customary to extend the domain of definitions to include 
the branch cuts, so that the function can be defined on the branch cut itself. There, 
the value of the function is taken as the limit of its values at points approaching the 
branch cut from one of its sides. In the numerical context, Kahan shows that using 
signed zeroes avoids having to make a choice [Kah87] . In the symbolic computation 
context a choice has to be made and the boundary of the domain is closed on one 
side and open on the other one. Making the choice of which side is the closed one 
is the adherence question [BBDP05]. For instance, the definition of In in [AS64] 
is taken so that ln(— 1) is i-K. Again, these choices have to be made consistently if 
several functions are involved. 

2.4. Inverse Functions. Inverse functions form a large class of functions that are 
commonly defined by analytic continuation. Suppose that / : C — > C is analytic, 
that f(xo) = 2/o and that f'(xo) ^ 0. Then there is a trivial function 

/ : {xq} -> {yo}, x y 

which clearly has an inverse. By the Inverse Function Theorem, this can be ex- 
tended to a neighbourhood of yo, and ultimately to the whole of C apart from those 
points where f'(x) = 0. 

Example 1. The basic example is defined as the inverse of 

f : C C, x^r x 2 . 
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Since /'(l) = 2 ^ 0, ice can define a function f in a neighbourhood of 1 
with = 1 and then extend it to larger connected domains. However, f 

cannot be continued arbitrarily far round the unit circle, for otherwise we would get 
the contradicting value / _1 (1) = — 1, see Section 2.1. 



Branch cuts for y/z, as well as other inverse functions like In z, z u , arcsin(2i), 
arccos(z), arctan(z), arcsinh(z), arccosh(z) and arctanh(z) are given in [AS64]. 
In all cases, they can be deduced from that of In once the function is expressed 
in terms of In (but even this expression is not neutral, as the second author was 
initially taught logarithms with a different branch cut). For these functions, Kahan 
claims [Kah87]: 

There can be no dispute about where to put the slits; their locations 
are deducible. However, Principal Values have too often been left 
ambiguous on the slits. 

In the terminology above, this means that the positioning question is soluble, and 
the problem is the adherence question. He states the following rules governing the 
location of the branch cuts: 

Rl. These functions / are extensions to C of a real elementary function analytic 
at every interior point of its domain, which is a segment S of the real axis. 

R2. Therefore, to preserve this analyticity (i.e. the convergence of the power 
series), the slits cannot intersect the interior of S. 

R3. Since the power series for / has real coefficients, f(z) = f(z) in a complex 
neighbourhood of the segment's interior, so this should extend throughout 
the range of definition. In particular, complex conjugation should map slits 
to themselves. 

R4. Similarly, the slits of an odd function should be invariant under reflection 

in the origin, i.e. z — > —z. 
R5. The slits must begin and end at singularities. 

While these rules are satisfied by the branch cuts of the inverse functions listed 
above, they do not completely specify their location, unless one adds a form of 
Occam's razor: 

R6. The slits might as well be straight lines. 
We shall interpret R4 in an extended way, by applying it as well when f(z)+f(—z) is 
a constant, as will be motivated by the example of inverse cotangent in Section 3.2. 

3.1. Worked example: arctan. Let us apply these rules to arctan, considered as 
the inverse of tan. Writing 



and solving a quadratic equation gives an expression for arctan in terms of In: 



From this expression one deduces that the singularities are located at ±i, so that it 
is analytic on R. Moreover, as the inverse of an odd function, arctan itself is odd. 
Hence we need a cut which 



3. Kahan's Rules 



arctan(z) = — - ( ln(l + ix) — ln(l — ix) 
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(R5): joins i and —i, 

(R3): is invariant under complex conjugation, and 
(R4): is invariant under z — > —z. 

(R6): We have the choice between a line from — i to i through and two lines 

—i — ti and i + ti, t > 0, meeting at infinity, but 
(R2): the first of these two options is not admissible, 
giving the classical branch cut z = + iy, \y\ > 1. 

3.2. The arccot dilemma. The strange case of arccot described in the introduction 
is still consistent with these rules. The key point is in Rl: in fact arccoti and 
arccotg were defined as different functions over R: they agreed on R + (in particular 
lim x ^, +ao arccoti (x) = lim^^+oo arccotg (a;) = 0), but not on R : 

arccoti(— 1) = 3-7r/4, 

arccotg(— 1) = — 7r/4. 

Therefore the limits at — oo arc different, and in fact arccotg is continuous at 
infinity (but discontinuous at 0). What should the branch cuts of these functions 
be? For arccoti, most of the reasoning of Section 3.1 applies. Strictly speaking, the 
function is not odd, but it is "odd apart from a constant", and hence the branch 
cuts should still be symmetric under z — > — z. Therefore it should have the same 
cuts as arctan, i.e. z — + iy, \y\ > 1. 

arccotg is odd, so all that reasoning applies, except that R2 no longer rules out 
the cut passing through 0. Indeed, since arccotg is discontinuous at 0, we are left 
with z = + iy, \y\ < 1 (the cut in [AS64, 9th printing]). 

4. Linear ordinary differential equations 

Many of the elementary, trigonometric, inverse trigonometric functions and hy- 
perbolic versions of those are part of the very large class of solutions of linear 
differential equations 1 (see Table 1). Here, we set to extend the previous set of 
rules to fix the location of the branch cuts in a way that is consistent with that of 
the previous section. Also, we only consider the case where the singularities are all 
regular singular points (meaning that the solutions have only algebraic-logarithmic 
behaviour in their neighbourhood). 

Let us assume throughout that we are given a linear ordinary differential equation 

n d l 

(3) L(y) = ^2,c i — i y = d, a,deC[x]. 

i=0 

At the cost of dividing by d, one differentiation and some re-normalisation, we can 
consider the homogeneous equivalent 

n+l 

(4) L( y ) = ^c i - 7 y = 0, Ci€C[x]. 

i=0 

More generally, we can homogenize Equation (3) whenever the inhomogeneous part 
d itself satisfies a linear o.d.e. of the form (4), as is the case with all examples in 
Table 1. Outside the zeros of c„ (or c n+ \), knowing y and sufficiently many of its 



Nonlinear equations have the major complication that it may not be obvious where the sin- 
gularities are, and indeed they may not be finite in number. Some entries of Table 1 do not satisfy 
a linear o.d.e., and this fact is indicated by a dash. 
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Table 1. Alternative definitions of functions 



Function 

exp 

log 

sin; cos 
tan; cot 
sec; esc 
arcsin; arccos 
arctan; arccot 
arcsec; arccsc 
sinh; cosh 
tanh; coth 
sech; csch 
arcsinh; arccosh 
arctanh; arccoth 
arcsech; arccsch 



Linear o.d.e. 

y' = y 
y' = i/x 
y" = -y 



v' = £1 



xyjx' 2 - 

y" = y 



y v^±T 

y = t^? 
v' = £1 



Definition by inverse 
log" 1 



arcsin 
arctan" 
arcsec 
sin" 
tan -1 
sec -1 
arcsinh -1 
arctanh - 
arcsech -1 
sinh -1 
tanh -1 
sech -1 



exp -1 

- ; arccos -1 

' -1 ; arccot -1 
1 



_1 ; arccsc 1 
• cos -1 
cot -1 

CSC -1 

arccosh - 1 
arccoth - 
arccsch -1 
cosh -1 
coth -1 
csch -1 



derivatives at some point xq (the obvious meaning of initial conditions) defines y 
as an analytic function in the neighbourhood of x - 

(5) y(x) = y(x ) + (x - x )y'{x ) H , 

where higher derivatives of y beyond the initial conditions are computed by applying 
(3) or (4) and their derivatives to the initial conditions. If it weren't for singularities, 
this would be an excellent definition. 

Example 2. The function can also be defined by 

xy'-^y = 0, y(l) = 1. 

Obviously, the leading coefficient x has a (regular) singularity at 0. 

4.1. Germs of branch cuts. In the vicinity of a regular singularity, the location 
of a branch cut can be shown by the form adopted for the local expansion. 

For instance, the branch cut for arctan joins i to — i along the imaginary axis via 
infinity (see Section 3.1). The local behaviour at i is therefore well described by 

—i 1 

(6) arctan(x) = — ln(l + ix) + iln\/2 + — (a? — *) + ■-■ 

written in such a way that the branch cut "heads north" . 

We can think of the precise formula used to encode the expansion at the singu- 
larity as encoding the germ of the branch cut, i.e. its local behaviour. The correct 
angle can always be achieved by rotating the argument. This solves the positioning 
problem as far as the germ of the branch cut is concerned. We also need to con- 
sider the adherence problem. Eq. (6) inherits the adherence from the logarithm, 
and therefore, for y > 1, means that 



arctan(0 + iy) = lim arctan(a; + iy) , 
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which is the adherence described in [Kah87] as "counter-clockwise continuity". 
When we need the other adherence, we simply use the fact that ln(l/x) = — ln(a;) 
except on the branch cut. 



4.2. Heuristic rules. The adaptation of Kahan's rules to an o.d.e. L(y) = 
together with a starting point is as follows: 

R2'. The branch cuts do not enter the circle of convergence. 

R3'. Complex conjugation is respected. 

R4'. Any symmetries inherent in the power series are respected. 

R5'. The branch cuts begin and end at singularities. 

R6'. The branch cuts are straight lines. 

R7'. The branch cuts are such that C less the branch cuts is simply connected. 

These subsume Kahan's rules, at the cost of explicitly requiring an initial value, 
which was implicit in his rules Rl and R2. He did not need an equivalent of R7' 
as his examples only had two singularities. In general, it is required so that the 
Monodromy Theorem (e.g. [Mar67, p. 269]) applies and guarantees uniqueness of 
function values. 

These rules do not necessarily completely determine the branch cut: a "random" 
differential equation with singularities scattered in the complex plane and no special 
symmetries will not be determined. Moreover, they do not give any guarantee of 
consistency between different functions. For instance, both functions g and q of (1) 
and (2) satisfy the same linear differential equations. Our rules that lead only to 
straight lines cannot be compatible with branch cuts that come from compositions 
of solutions of simpler differential equations with algebraic functions. However, 
they serve the simple purpose of producing useful and correct branch cuts in a wide 
variety of cases, including all those discussed before. 



4.3. Worked example: arctan. We apply these rules to arctan, now defined by 

The singularities of this differential equation are clearly at x = ±i, and the function 
so defined is odd. Hence we need a cut which: 

(R5'): joins i and —i, 

(R6'): does it in a straight line, 

(R3'): is invariant under complex conjugation, 

(R4'): is invariant under z — > — z, 

(R2'): does not enter the unit disk. 

Thus we find again the classical branch cut z = + iy, \y\ > 1. We then deduce 
expansions at the singularities that match the germs of this cut as in (6). 



4.4. 5 In a; or lncc 5 ? Once In has been defined, the functions F\ — 51nx and 
F 2 = In a; 5 are different: Fi(i) = ^ while F 2 (i) — ^. Nevertheless, they are 
both solutions to xy' — 5 = with y(l) = 0. Our approach would make the 
choice 51nx with only one branch cut, while F 2 has five branch cuts, at angles of 
{1, 3, 5, 7, 9}7r/10, thus making the domain not connected (and violating R7'). 
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4.5. A harder example. Let us consider the functions / defined by 



J l + x 4 ' 
or, if one prefers homogeneous equations, 

x{l + x 4 )f" + (3a; 4 - f )/' = 0. 



In both cases, we assume we are given real initial conditions at 0. This ex- 
ample is selected because it has two simple, linearly independent solutions — 1 
and arctan(x 2 ) — to compare with the result, but the method does not use this 
information and would apply even if no such solution could be found. So here is 
what we get: 

(R5'): The equation has four regular singularities at 



(one can check that is just an apparent singularity by exhibiting a basis of 
formal power series solutions and that oo is not a singularity by changing x 
into 1/x). 

(R6'): These four singularities have to be connected by straight lines. 
(R2'): We cannot connect the singularities pairwise (in either way!) without 
going to infinity. 

(R4'): The symmetry f(ix) = —f(x) can be checked directly from the equa- 
tion, so that branch cuts should be mapped to branch cuts by a rotation 
of tt/2. 

(R3'): Reality implies that branch cuts are also mapped to branch cuts by 
horizontal symmetry. 
We are thus left with only the following choice: Cuts that "head northeast" from 
"northwest" from etc., all meeting at infinity. This is indeed consistent 

with arctan(x 2 ). 

It is worth noting that this function actually also admits branch cuts that violate 
R2' and R7': for example we can connect to ± ^ L , and =^j± to ± ^-- This 

is a peculiarity of our construction, and the fact that these are valid follows, not 
from the Monodromy Theorem, but from the fact that the residues at these branch 
points are equal and opposite. 



When it comes to converting an analytic (be it linear ordinary differential equa- 
tion, inverse function, or possibly other) definition of a function into a well-defined 
single-valued one, so that one can answer questions such as "what is ln(— 1)?" or 
"what is arctan(2i)" , branch cuts may need to be imposed on the locally analytic 
function. While the definition of the function may stipulate the endpoints of the 
cut, it docs not, in general, specify the location of the cut between its endpoints, 
nor indeed even the germ of the cut at the singularities. 

We have given a simple set of rules that is convenient when nothing else is 
known about the function. This set of rules is sufficient to recover the classical 
branch cuts of the elementary inverse trigonometric or hyperbolic trigonometric 
functions. However, it is important to remember that this is only a useful heuristic, 
while there are cases where different cuts are dictated by the application. In a 
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specific context, getting the right overall function defined by a formula can even 
require inconsistent choices of the branch cuts of component functions: see, e.g., 
the Joukowski map studied in [Hen74, pp. 294-8]) and reported on in [DavlO]. 
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